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Рассмотрены проблемы, методы и результаты создания программ имитационного моделирования, а также их 
алгоритмического и программного обеспечения. Проблемы рассмотрены на примерах технологических рас- 
пределённых объектов, таких как участок магистрального газопровода и многослойная структура материа- 
лов, в которой протекают тепловые процессы. 
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Введение. Решая многие задачи оптимизации, диагностики, управления технологическими про- 
цессами, необходимо получать подробную информацию об этом процессе. Как правило, наиболее 
важная для этих задач информация носит «предсказывающий» характер. Такая информация по- 
рождается в ходе имитации, она позволяет ответить на вопрос, как вёл бы себя объект управле- 
ния и какие бы там протекали процессы, если определённым образом воздействовать на объект. 
Обладая такой информацией, можно строить эффективные и близкие к оптимальным технические 
системы. 

Имитация в технических системах находит применение в информационно-управляющих 

диалоговых системах, которые сигнализируют и блокируют последствия «ошибочных» действий 
оператора технологического процесса. 
Имитационное моделирование в информационной поддержке технологических про- 
цессов. Особую роль имитация технологических процессов приобретает, когда поведение управ- 
ляемой системы чрезвычайно сложно предсказать. Одной из таких систем является система рас- 
пределенных технологических объектов (РТО). Современные специализированные программные 
пакеты справляются с анализом и синтезом систем РТО на основе их статических характеристик. 
Непосредственно РТО — широкое понятие, описывающее технологические агрегаты, имеющие де- 
ло с распределенной в пространстве средой, в то же время под РТО иногда понимают комплекс 
технологических объектов, пространственно распределённых по территории производства. В 
дальнейшем примем первый смысл этого выражения. Наиболее сложные математические модели 
получаются для распределённых технологических объектов, в которых процессы протекают не 
только во времени, но и в пространстве. К ним относятся: трубчатые и колонные реакторы, печи, 
накопители газа (газгольдеры), газопроводы, транспортеры, механические фильтры и фильтруе- 
мый материал, среды передачи информации электромагнитными волнами. Распределённость учи- 
тывается в математических моделях технологических объектов, когда объем, занимаемый рабо- 
чей средой, достаточно большой по сравнению с величинами локальных неоднородностей (волн) 
в переходных процессах. 

Имитационное моделирование как надстройка над математическим моделированием при- 
меняется в тех случаях, когда методы аналитического вывода решения оказываются ограничены 
выбором пространств и метрик, а также необходимой строгостью выводов и доказательств. Суть 
имитационного моделирования заключается в образовании таких информационных связей 
(рис. 1) в уравнениях, которые бы имитировали связи в моделируемом объекте другой природы, 
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например энергетические, массовые и др. Помимо связей в имитационных моделях информация 
обрабатывается в узлах по более сложным законам, задаваемым уравнениями, иногда дифферен- 
циальными. Наиболее просто имитационная модель реализуется в виде программы на языке вы- 
сокого уровня или как структурно-блочная схема в среде специализированных пакетов (Апу5У$, 
Ном/\Ляюп, СОМФО(, ЗитийтК, и т. д.). 
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Рис. 1. Имитационная модель газоперекачивающего агрегата 
в среде пакета 5ититпк 


Проблемы имитационного моделирования. В имитационном моделировании применяются 
системы методов численного решения, применение которых, в свою очередь, связано с приняти- 
ем системы обоснованных допущений, позволяющих исследование объекта и процессов, там про- 
текающих, свести к исследованию функций и числовых рядов. Вариации допущений в строну уси- 
ления дают простые математические модели, позволяющие предсказать поведение объекта за 
небольшое число вычислительных операций. Ослабляя допущения, в ходе аналитического моде- 
лирования исследователь-разработчик получает более точную модель. Такие модели описывают- 
ся более сложными системами уравнений, требующими большего количества вычислений. Пред- 
сказывающая способность таких формул по отношению к процессу будет лучшей. Поэтому на эта- 
пе построения вычислительных формул и самой имитационной модели стоит проблема принятия 
системы допущений, которая, с одной стороны, позволит получить нужную информацию, а с дру- 
гой стороны, позволит рационально организовать вычислительный процесс. 

При моделировании многих распределённых технологических объектов принимают допу- 
щения, которые условно можно разделить на необходимые и дополнительные. К необходимым 
можно отнести параметры функций состояния (интенсивности их изменения известны только на 
некоторых границах или сечениях); среда в технологическом аппарате локально ведёт себя во 
всех точках достаточно большой области одинаково, а погрешности математических аппроксима- 
ций и пространственно-временной дискретизации качественно не изменяют поведение модели по 
отношению к моделируемому объекту. Эти допущения позволяют выписать математическую мо- 
дель. Допущения, принимаемые для упрощения математической и имитационной модели и, как 
следствие, для уменьшения объёма вычислений, можно отнести к дополнительным. Описанные 
зависимости поясняются далее. 

Желание исследователя сделать моделирование максимально точным приводит к умень- 
шению глубины допущений, принимаемых при математическом моделировании, что приводит к 
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появлению сложностей математической и алгоритмической природы. Большое количество допу- 
щений позволяет получать простые связи между причиной и следствием выраженными в виде ал- 
гебраических уравнений; при меньшем числе допущений получают математические модели в ви- 
де системы дифференциальных уравнений. Невозможность в большинстве случаев принять до- 
пущение о сосредоточенности процессов в распределённых объектах приводит к математическим 
моделям в виде систем дифференциальных уравнений в частных производных; при ещё меньшей 
глубине допущений можно получить интегродифференциальные и статистические уравнения. Ес- 
тественно, что с первыми объектами — алгебраическими уравнениями — проще иметь дело, а 
дифференциальные уравнения в частных производных вызывают определённые математические 
и алгоритмические сложности, особенно когда некоторые математические зависимости, связы- 
вающие параметры среды -— газа, жидкости — являются нелинейными. Ярким примером может 
служить модель участка магистрального газопровода (УМГП), в которой учёт распределённости 
процессов в пространстве и времени — естественная необходимость, а допущения, которые 
рационально принять: течение осесимметричное, допустимо радиальное усреденение параметров 
при турбулентном движении (рис. 2) [1]. В этом случае в качестве модели выступает «компакт- 
ный» математический объект — система дифференциальных уравнений в частных производных 
[2, 3]: 
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В (1) Е х- независимые аргументы дифференциального уравнения -— время и пространственная 
координата; р, у, Т — неизвестные функции двух координат: давление, скорость и температура; 
и, К, Г- параметры транспортируемого газа: молярная масса, универсальная газовая постоянная, 
количество степеней свободы молекулы газа; ^, О — параметры, характеризующие трубопровод: 
коэффициент трения (параметр Дарси — Вейсбаха), диаметр поперечного сечения газопровода; 
Т‹р - характеристика среды залегания трубопровода — температура, характеризующая интенсив- 
ность потери тепла и, как следствие этого, — полной энергии газа. 





Рис. 2. Расчетная схема взаимодействия частиц газа 


Дальнейшие допущения ещё больше огрубляют математическую модель. Предположение 
об изотермическом поведении газа при транспортировке даёт модель (2), а исключение из моде- 
ли динамических процессов даёт модель статики (3): 
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На выборе системы допущений сложности и проблемы имитационного моделирования не 
заканчиваются. Необходимо исследование математических моделей и решение уравнений, их 
представляющих. Количество итераций вычислительного алгоритма, как правило, велико и зави- 
сит от требуемой точности и от пространственно-временных масштабов моделируемого объекта. 
Данную проблему осложняет ещё одно «необычное» свойство большинства исследуемых процес- 
сов — жёсткость. Дифференциальные уравнения, особенно учитывающие несколько видов взаи- 
модействий (например, химические и механические), часто имеют жёсткие решения [4]. Жёст- 
кость — возможность возникновения в модели взаимодействующих процессов с существенно раз- 
ными скоростями протекания. Область решения более простых алгебраических уравнений может 
быть «сильно овражистой», что может потребовать большого количества итераций для численно- 
го решения. Жёсткость математической модели транспортировки газа по участку магистрального 
газопровода иллюстрируется спектральными диаграммами линеаризованного приближения мате- 
матической модели (2). На графиках рис. 3 видно расслоение мод на быстрые и медленные. 

Неоднозначен и слабо формализован подход к дискретизации распределённых объектов 
при аппроксимации дифференциальных уравнений в частных производных (ДУЧП). Возможен 
подход с полной аппроксимацией. В результате получается система алгебраических уравнений 
большой размерности, описываемая разреженными матрицами. Желание авторов программ мак- 
симально использовать хорошо зарекомендовавшие себя разностные методы приводит к частич- 
ной аппроксимации ДУЧП по одной из координат, например пространственной, в результате по- 
лучается система дифференциальных уравнений большой размерности, решаемая в дальнейшем 
такими методами, как метод Рунге - Кутты, Адамса — Крылова, Милна. Такой подход называется 
методом прямых или в силу поэтапного перехода — полудискретизацией. 
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Рис. 3. Расположение мод линеаризованной модели на комплексной плоскости. Цифрами показаны моды, 
полученные дискретизацией по пространству на соответствующее количество элементов - ячеек 


794 


Вестник ДГТУ. 2011. Т. 11, № 6(57) 


Для вывода, построения алгоритмов моделирования строятся вычислительные формулы, 
основанные на той или иной математической модели. Для вывода вычислительных формул ис- 
пользуется прямая ячеечная алгебраическая аппроксимация дифференциальных уравнений. Она 
применяется в методах конечных разностей, конечных объёмов (рис. 4) [5-7], конечных элемен- 
тов и в подавляющем большинстве случаев нарушает интегральный баланс массы, энергии и им- 
пульса. Нарушение баланса можно компенсировать построением неравномерной сетки (рис. 5). 
Этот вопрос подробно исследовался в кандидатской диссертации автора [8] и других публикациях 
[9] о модели участка магистрального газопровода. Для уравнения неразрывности (4) был получен 
закон формирования пространственной сетки (5): 
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Рис. 4. Схема полудискретизации пространства для моделирования участка магистрального газопровода 
методом конечных объёмов 
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Рис. 5. Диаграммы распределения импульса 
при дискретизации параметров в уравнении неразрывности 


Сохранение подобного баланса называют консервативностью [6]. Одним из способов по- 
лучения консервативных разностных схем и моделей является непосредственная аппроксимация 
(по квадратурным формулам) интегралов энергий исходной математической модели. 
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Для повышения точности моделирования используют разностные схемы повышенного по- 
рядка (второго и выше) точности. Такие схемы неплохо моделируют медленноменяющиеся коле- 
бания и «гладкие» зависимости времени. При попытке воспроизвести резкие изменения в модели 
такие разностные схемы ведут себя излишне колебательно, порождая паразитные гармоники. Для 
решения этой проблемы возвращаются к разностным схемам первого порядка: правосторонним и 
левосторонним, а общее решение этой проблемы приводит к разностным схемам со сложным пе- 
реключательным шаблоном, называемым схемой с уменьшением полной вариации (ТУБ — То! 
\апаноптп Оиттйюоп). Выражения, которые их описывают, называются лимитерами. В современных 
научных работах по моделированию течений различной природы исследовано множество лими- 
теров (ИтКег): «тт-по@», «уап [еег», «зирегрее» и др. [10]. 

Замена разностного шаблона пространственной аппроксимации с правостороннего на ле- 

восторонний шаблон первого порядка (две точки) сильно сказывается на решениях обыкновенно- 
го уравнения переноса. Для этого уравнения противопотоковые дают устойчивое решение [6], 
что адекватно физике процесса. При выполнении шагов по времени в разностных выражениях 
соблюдается такая же зависимость. Неявные разностные схемы абсолютно устойчивы при любых 
шагах численного метода. Явные разностные схемы устойчивы условно, т. е. величина шага огра- 
ничена условием Куранта — Фридрихса — Леви (КФЛ), что обосновывается спектральным анали- 
зом, при этом условие КФЛ называется спектральным условием Неймана [6, 7]. В кандидатской 
диссертации [7] подробно изучены спектры линеаризованного приближения математической мо- 
дели магистрального газопровода, которые показывают предпочтительность центрально- 
разностной схемы аппроксимации и левосторонней схемы. 
Программное обеспечение, поддерживающее исследования. Разработано программное 
обеспечение (ПО), осуществляющее моделирование участка магистрального газопровода с парой 
компрессорных станций. ПО представлено в виде законченного модуля — приложения ОС \М/паом5 
и программы для среды пакета МаЧаб, которая в силу удобного синтаксиса системы МаЧаБ ис- 
пользуется как экспериментальный образец, в котором легко перестраиваются информационные 
связи, в связи с чем удобно получить модель кольцевого газопровода, энергия газа в котором 
поддерживается одной компрессорной станцией. Вычислительные выражения получены нами ме- 
тодом полудискретизации исходной математической модели [11]. 

Для сравнительного анализа подходов к моделированию была разработана программа мо- 
делирования одномерного температурного поля набора слоёв материалов методом элементарных 
энергетических балансов (рис. 6). Программа обладает графическим интерфейсом и работает под 
управлением системы МаНаБ, используя среду СУТВЕ [12, 13]. 
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Рис. 6. Интерфейс программы моделирования одномерного температурного поля 
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В ходе исследования проблем и походов к моделированию распределённых объектов была 
поставлена задача подробного исследования отдельных этапов моделирования. В соавторстве со 
студентом Д. Лауром, выполняющим дипломный проект, был разработан программный комплекс 
исследования разностных аппроксимаций распределённой модели участка газопровода (рис. 7). 
Допускается ввод произвольных разностных шаблонов и выполнение моделирования по времен- 
ным шагам. 
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Рис. 7. Интерфейс интегрированной среды исследования разностных аппроксимаций 


Отдельным направлением ведётся разработка программы имитационного моделирования 
типовых распределённых объектов, которая нужна для того, чтобы верифицировать алгоритмы 
моделирования и проверить адекватность моделей технологическим объектам [14]. 

В ходе разработки были выработаны принципы построения программ имитационного мо- 

делирования распределённых технологических объектов общей природы. Были выделены слои 
описания моделей, опубликована иерархия классов и архитектура взаимодействия объектов в хо- 
де имитационного моделирования [8, 15]. 
Заключение. В рамках исследования построены алгоритмы, позволяющие формировать про- 
странственную сетку сбалансированной аппроксимации для шаблонов разностных схем с произ- 
вольной структурой; разработаны и протестированы алгоритмы численной интерполяции, интег- 
рирования и дифференцирования, позволяющие моделировать транспортировочные процессы на 
разнесенных сетках произвольными разностными шаблонами. Разработанные алгоритмы позво- 
ляют имитационными методами подобрать рациональный по точности и сходимости разностный 
аппроксимирующий шаблон. 
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РВОВЕЕМ$ АМО 5ТМОГАТТОМ МЕТНОО$ ОЕ ОТЪЗТКТВОТЕО ОВЗЕСТЗ 


М.У. КООТМОУ 
(Роп Чае Тесйптса! Утмег$Ку) 


ргоБетх, тетоа$ апа гезиЁ5 оЁ те ти/айоп ргодгатз деуеортейе а$ и/е! аз {пе а/догтис зирро!Е апа 
5ОНуаге аге сопз!аегеа. Тпе ргоМет$ аге @5сиз5еа т {епт о! {пе есбппоод/с @триеа обес, 5иср аз те 
тат даз р/ре/те а/5оп апа {пе таёепа! ти!ауегеа $гисйиге ипеге {пегта! ргосез$е5 ке р/асе. 

Кеуигога$: Оотршеа оБлесЕ5, эти/айоп (есптие, татетайса! этиайоп, хтиайоп ргодгатз. 
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